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Abstract. In this survey article, we investigate the spectral properties of 
fractal differential operators on self-similar fractals. In particular, we discuss 
the decimation method, which introduces a rcnormalization map whose dy- 
namics describes the spectrum of the operator. In the case of the bounded 
Sicrpinski gasket, the renormalization map is a polynomial of one variable on 
the complex plane. The decimation method has been generalized by C. Sabot 
to other fractals with blow-ups and the resulting associated renormalization 
map is then a multi-variable rational function on a complex projective space. 
Furthermore, the dynamics associated with the iteration of the renormalization 
map plays a key role in obtaining a suitable factorization of the spectral zeta 
function of fractal differential operators. In this context, we discuss the works 
of A. Teplyaev and of the authors regarding the examples of the bounded and 
unbounded Sicrpinski gaskets as well as of fractal Sturm-Liouville differential 
operators on the half-line. 
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1. Introduction 



From the probabilistic point of view, the Laplacian on the Sierpinski gasket SG 
was introduced independently by S. Goldstein in |18| and S. Kusuoka in |27j (and 
a little later, by M. Barlow and E. Perkins in [6]), as the generator of the semigroup 
associated with Brownian motion on SG. (See, e.g., [4j and [5] for early reviews of 
the subject of diffusions and random walks on self-similar fractals.) However, from 
the point of view of analysis, which will be our main concern here, the Laplace 
operator was first defined by J. Kigami |23j for the Sierpinski gasket and was later 
extended in |24j to a class of self-similar fractals, called the post critically finite 
sets (p.c.f. sets). (See, e.g., (25j and |51j for a detailed exposition.) The Laplacian 
on a p.c.f. set is defined as the limit of a sequence of Laplacians of finite graphs 
that approximate the fractal. Following the work of the physicists R. Rammal and 
G. Toulouse @0j[4T], M. Fukushima and T. Shima [THUD] studied the eigenvalue 
problem associated with the Laplacian on the Sierpinski gasket and introduced the 
decimation method in order to give an explicit construction of the set of eigenvalues. 
The decimation method, described in §2 of the present paper, is a process through 
which we find the spectrum of the Laplacian on a fractal (in a certain class of self- 
similar sets) via the iteration of a rational function of a single complex variable, 
called the renormalization map. In the case of the finite (or bounded) Sierpinski 
gasket, this rational map is a polynomial on the complex plane and its dynamics is 
not too difficult to understand in order to analyze the spectrum of the Laplacian. 

Later on, C. Sabot (|43j-[47j) generalized the decimation method to Laplacians 
defined on a class of finitely-ramified self-similar sets with blow-ups. First, in 
|43L I44L I47| . he studied fractal Sturm-Liouville operators on the half-line, viewed 
as a blow-up of the self-similar unit interval, and discovered that the corresponding 
decimation method then involves the dynamics of a rational map which is no longer 
a function of a single complex variable but is instead defined on the two-dimensional 
complex projective space; see §3. (It therefore arises from a homogeneous rational 
function of three complex variables.) This rational map is initially defined on a 
space of quadratic forms associated with the fractal and its construction involves 
the notion of trace of a symmetric matrix on a (finite) subset of this space; see §4. 
From the point of view of multi-variable complex dynamics, the set of eigenvalues 
(or spectrum) is best understood in terms of an invariant curve under the iteration 
of the rational map. 
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Using as a model the unbounded (or infinite) Sierpinski gasket (the so-called 
Sierpinski lattice, see §4 and Figure[4]), Sabot |46j used Grassmann algebras in order 
to construct the renormalization map for other lattices based on (symmetric) finitely 
ramified self-similar sets. The idea is that one embeds the space of symmetric 
matrices in a Grassmann algebra in order to analyze the operation of taking the 
trace on a (suitable) finite subset. In some sense, this enables one to linearize 
this operation. The polynomial associated with the classical bounded Sierpinski 
gasket, initially introduced in the work of Rammal and Toulouse |40L 141] and later 
rigorously formalized in |17L 150] , can then be recovered from the renormalization 
map associated with the unbounded Sierpisnki gasket. 

Finally, in recent work, Lai and Lapidus [28j have studied the spectral zeta 
function of the Laplacian on a suitable self-similar set and established a factorization 
formula for the associated spectral zeta function in terms of a certain hyperfunction, 
a geometric zeta function and a zeta function associated with the iteration of a 
renormalization map, which is a multi-variable rational map acting in a complex 
projective space. This latter work in |28j extends to several complex variables an 
earlier factorization formula due to Teplyaev |54|, 155] , itself extending the second 
author's factorization formula |30l 13 lj for the spectral zeta function of a fractal 
string (see also |33l 1341 136] for various applications of, and motivations for, the 
latter factorization). We survey some of these results in the last part of this paper; 
see §5. 

In closing this introduction, we mention that the work of |28j described in §5 
focuses on two different models, namely, fractal Sturm-Liouville differential oper- 
ators on the half-line (as in §3) and the infinite (or unbounded) Sierpinski gasket 
(as in §4). In each of these cases, the Dirac delta hyperfunction plays a key role in 
the rigorous mathematical formulation of the factorization results. 

2. The bounded Sierpinski gasket 

The (bounded or finite) Sierpinski gasket (SG) is generated by the iterated 
function system (IFS) consisting of three contraction mappings <F,- : M 2 — > W 2 
defined by 

(2.1) $ j (x) = ^{x-q j ) + q j 

for j — 0, 1, 2, where qo, q\, q2 are the vertices of an equilateral triangle. (See Figure 
1.) Note that each $j has a unique fixed point, namely, qj. The Sierpinski gasket is 
the unique (nonempty) compact subset of K 2 such that SG = <&i(SG) U <fr 2 (SG) U 
MSG). 

For each integer m > 0, denote by T m the mth level finite graph approximating 
SG, and let r be the complete graph on Vq :— {q\, q 2 , 93}- The set V m of vertices 
of r m is defined recursively as follows: 

2 

Vm ■= \J$j{V m -i), 
3=0 

for each integer m > 1. 

Let £ 2 (V m ) be the space of real- valued functions on V m , equipped with the 
standard inner product (it, v) — J2 x ev u{x)v(x). The discrete Laplacian on £ 2 {V m ) 
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Figure 1. The (bounded) Sierpinski gasket SG. 

(or the finite graph Laplacian on L m ) is defined by 

A m /(aO =*£/(*)-/(</), 

where x G T m and the sum is extended over all neighbors y of x in the graph T m . 

We then define the Laplacian A = A M on SG as the following limit of rescaled 
finite- difference operators: 

A f{x) = lim 5 m A m f(x). 

m— >oo 

The factor 5 is the product of the scaling factor 3 for the natural Haursdorff measure 
on SG and the renormalized factor | for the energy. Here and thereafter, when we 
write A = A M , the subscript /j, refers to the natural self-similar (or equivalently, in 
this case, Hausdorff) probability measure on SG. 

The graph energy on each V m (to > 0) is defined by 

£ m (u,u)=(^j ^2(u(x) - u(y)) 2 , 

which docs not change under the process of harmonic extension, and the graph 
energy on SG is then defined by 

£(u, u) = sup £ m {u, u) = lim £ m (u,u) 

m>0 m-HX 

(this limit always exists in [0, oo] since the sequence {£ m (u, m)}^ =0 i s nondecreas- 
ing). In the sequel, we also write £(u) — £(u, u), in short, to refer to this quadratic 
form. (By definition, u belongs to the domain of £ if and only if £{u) < oo.) The 
associated bilinear form £(u,v) is then defined by polarization. 

Suppose u is a function defined on Vo, with values at each of the three vertices 
of Vo denoted by a, b, and c; see Figure 2. We want to extend u to V\ in such a 
way that it minimizes the energy. 

Let u be the harmonic extension of u to V\ and denote by x, y and z the values 
of u to be determined at each of the three vertices of V\ \ V ; see Figure 2. Since, 
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Figure 2. The values of the harmonic extension u at each vertex 
in V\. The values a, b, c at each vertex in Vo are prescribed, whereas 
the values x, y, z at each vertex in V\ \ Vq are uniquely determined 
by requiring that the energy be minimized in the passage from Vo 
to Vi. 



by definition, u minimizes £\{v) subject to the constraint v = u on Vo, we can take 
the partial derivatives with respect to x, y, z and set them equal to zero to obtain 

Ax = b + c + y + z 
4y = a + c + x + z 
Az = a + b + x + y. 

These equations express the mean value property of a (discrete) harmonic func- 
tion, according to which the function value at each of the junction points is the 
average of the function values of the four neighboring points in the graph. We can 
use the matrix representation of these equations 




in order to obtain the following solutions: 



1 

x = -a 
5 



y=- 5 a 



z = -a - 
5 



2 , 2 

1 r. 2 

7 b+ 7 c 
5 5 

2, 1 



The harmonic extension u therefore satisfies the \ — % rule. See 

5 5 

|51j . (More generally, this rule also holds for the harmonic extension on each m- 
cell of SG.) More specifically, iterating this process to each finite graph V m , one 
then defines u on the countable set of vertices V* :— ii m >oV m and finally, extends 
it to all of SG, by continuity and using the density of V* in SG. The resulting 
function, still denoted by u, is called the harmonic extension of u. According to 
Definition 2.1 below (which is also a theorem), it is the unique harmonic function 
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(i.e., Am = 0) such that u = u on Vq. Equivalently, it is the unique minimizer of 
the energy functional £ = £(v) subject to the constraint v = u on Vq. 

Definition 2.1. A harmonic function on SG {with boundary value u on Vq) is 
a continuous function whose restriction to any T m or V m is the harmonic extension 
of u. In other words, it is the unique solution of the following Poisson problem: 
Av = 0, v = u on Vq. It must therefore necessarily coincide with the harmonic 
extension u of u to all of SG. 

2.1. Spectral properties of the Laplacian on the Sierpinski gasket. 

The Laplacian operators on p.c.f. self-similar fractals are defined similarly via a 
suitable approximation. To study the spectrum of the Laplacian, we consider the 
equation — Am = Am, where u is a continuous function. The spectrum of the Lapla- 
cian on the Sierpinski gasket was first studied in detail by the physicists R. Rammal 
and G. Toulouse [40l |41] . Later on, M. Fukushima and T. Shima [T71 [50] gave 
a precise mathematical description of the eigenvalues and the eigenfunctions. Still 
in the case of the Sierpinski gasket, Rammal and Toulouse discovered interesting 
relations between the spectrum of the discrete Laplace operator and the iteration 
of a polynomial of one complex variable, R — R{z) := z(5 — 4z). More precisely, 
for any m > 0, if A is an eigenvalue of — A m+ i on r m+ i, then A(5 — 4A) is an 
eigenvalue of — A m on L m . Thus, the relationship between the eigenvalues of the 
Laplacians on one graph and it successor can be described by a quadratic equation, 
A m = A„ l+ i(5 — 4A m+ i) = i?(A m+ i). The restriction to V m of any eigenfunction 
belonging to A m+ i is an eigenfunction belonging to A m . The relationship between 
the eigenvalues A m and A„ l+1 of — A m and — A m+1 , respectively, can be found by 
comparing the corresponding eigenvalue problem for a point common to both V m 
and V m+ i. 

Theorem 2.2 (Fukushima and Shima, [17], [50] ). 

(i) If u is an eigenfunction of — A m+ i with eigenvalue A (that is, — A m+ iM = Am), 
and if A ^ B, then — A m (M|y m ) = R(X)u\v m , where B :— {\, |, |} is the set 
of 'forbidden' eigenvalues and u\v m is the restriction of u to V m . 
(ii) If — A rn u — R(X)u and A ^ B, then there exists a unique extension w of u to 
V m +i such that —A m+ iw = Xw. 

At any given level m, there are two kind of eigenvalues of — A m , called the initial 
and continued eigenvalues. The continued eigenvalues arise from the spectrum of 



— A m _i via the decimation method (described in Theorem 2.2), and the remaining 
eigenvalues are called the initial eigenvalues. The forbidden eigenvalues {\, |, |} in 



Theorem 2.2 have no predecessor, i.e., they are the initial eigenvalues. Furthermore, 
the exclusion of the eigenvalue \ can be explained by showing that \ is an eigenvalue 
of — A m only for m = 1. (See Figure [3j) 

Given that the eigenvalues of — A are {0, |}, we consider the inverse images 

of and | under R (that is, their images under R-(z) — 5-v/2 g~ 16 - and R+(z) = 

5+v/2 |~ 16 -, the two inverse branches of the quadratic polynomial R(z) = z(5 — 4z)), 
to obtain the eigenvalues of — Ai. The continuation of this process generates the 
entire set of eigenvalues for each level. The diagram provided in Figure [3] illustrates 
the eigenvalues associated with each graph Laplacian — A m , for m = 0, 1,2, in 
terms of the inverse iterates of the polynomial R. 




Figure 3. The eigenvalues of SG are obtained via the decimation 
method as limits of inverse images of the renormalization map R. 
At each level, the two branches of the tree correspond to the two 
inverse branches of the quadratic polynomial R, denoted by R + 
and R- . The number in the bracket represents the multiplicity of 
the corresponding eigenvalue. 



The spectrum of the Sierpinski gasket is the renormalized limit of the spectra 
of the graph Laplacians — A m . More specifically, each eigenvalue satisfying the 
equation — A^u — Xu can be written as 



(2.2) 



A = 



lim 5* 

n—>ao 



for a sequence {X m }™ =mo such that X m = A m+1 (5 - 4A. m+ i) = R(X m ) for all 
m > tuq and for some smallest integer mo (which is allowed to depend on A). Note 
that for m > tuq, X rn does not coincide with any of the forbidden eigenvalues in 
{|, |, |}, whereas X mo belongs to the set {|, |, |}. Furthermore, the values A m are 
determined by the solutions of A m = A m+ i(5 — 4A m+ i): A m +i = 5 + £ '"V / y~ 16A "' ; 
where e m = ±1, provided that the limit in ( |2.2| exists. The limit A only exists 
if e m = —1 for all but finitely many integers m. In that case, A is an eigenvalue 
of — A M where, as before, A = A M denotes the Laplacian on SG. And conversely, 
every eigenvalue of — A M can be obtained in this manner. 

In the next section, we will discuss a generalization (due to Sabot \44\ 145] ) 
of the decimation method to rational functions of several complex variables, when 
presenting the case of fractal Sturm-Liouville differential operators. The extended 
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(multi-variable) decimation method is not valid for all self-similar fractals, in gen- 
eral. However, it does apply to a large class of symmetric p.c.f. self-similar fractals 
(whereas the original single- variable decimation method only applied to a rather re- 
stricted and difficult to characterize class of symmetric finitely ramified and hence, 
p.c.f., self-similar fractals). Therefore, this generalization is a very significant exten- 
sion of the original (single-variable) decimation method, for which it also provides 
a nice geometric and algebraic interpretation (see §4). 

In fact, mathematically, the truly beautiful underlying structure of the extended 
decimation method is only revealed by considering the multi-variable case, even in 
the original setting of a single complex variable. (The latter one-variable case 
should really be formulated in terms of two complex variables or equivalently, on 
the complex projective line, and via phase space symplectic geometry or Grassmann 
algebras, in terms of two conjugate variables. In hindsight, the reduction to a single 
variable is simply a confusing, albeit convenient, artifact in this situation.) 

Remark 2.3. The interested reader can find in O [3] some detailed and rel- 
atively elementary computations pertaining to the decimation method and the as- 
sociated renormalization map (viewed only as a rational function of one complex 
variable) in the case of certain examples of finitely ramified, symmetric self-similar 
fractals. 

3. Generalization of the decimation method 

3.1. The fractal Sturm Liouville operator. C. Sabot, in a series of papers 
( 43 1 [47 J, extended the decimation method to Laplacians defined on a class of 
(symmetric) finitely ramified (really, p.c.f.) self-similar sets with blow-ups. This 
extension involves the dynamics of rational functions of several complex variables. 
We discuss the prototypical example he studied, fractal Laplacians on the blow-up 
I<oo> — [0, oo) of the unit interval / = /<o> = [0, 1]. From now on, we will assume 
that 

a 1 

(3.1) < a < 1, 6=1 — a, 5 = , and 7 = — 

1 — a a(l — a) 

Consider the contraction mappings from / = [0, 1] to itself given by 

*i(x) = ax, ^ 2 (x) = 1 - (1 - a)(l - x), 

and the unique self-similar measure m on [0, 1] such that for all / G C([0, 1]) (the 
space of continous functions on / = [0, 1]), 

(3.2) f fdm = b [ f o * x dm + (1 - b) [ foV 2 dm. 
Jo Jo Jo 

Here, I is viewed as the unique self-similar set (in the sense of |21| ) associated with 
the iterated function system {^1, ^2}- 



(3.3) I = *i(/)U* a (-0. 

Define H <0> = — r-jr, the free Hamiltonian with Dirichlet boundary condi- 
tions on [0, 1], by H < Q > f = g on the domain 

{/ 6 L 2 (I, m) : 3g E L 2 (I, m),f(x) = cx+d+ J jT g(z)dm(z)dy, /(0) = /(l) = oj 
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The operator i?<o> is the infinitesimal generator associated with the Dirichlet form 
(a, T>) given by 

a(f,g)= f f'g'dx, for 

Jo 

where 

V = {f eL 2 (I,m): f EL 2 (I,dx)}. 

As can be easily checked, the Dirichlet form a satisfies the similarity equation 

(3.4) a(f) = a - x a{f o x ) + (1 - a^aif o # 2 ), 

where we denote the quadratic form a(f,f) by a(/). (See, e.g., [13] for an exposi- 
tion.) 

The idea is that viewing the unit interval I as a self-similar set, as in (3.3) 
above, we construct an increasing sequence of intervals I<n>, f°r n = 0, 1, 2, by 
blowing-up the initial unit interval by the scaling ratio a~ n . Hence, we can extend 
all the objects involved, (m, a, H <0> ), from I = I<o> to /<„> = \&T n (J) = [0, a - ™], 
which can be expressed as a self-similar set as follows: 

I<n> = IJ **!...*„(/<«>), 

where i n ) G {1, 2}™. Here, we have set ^ x ...i n = ° ••• ° 

More precisely, for each n > 0, we define the self-similar 

7<«> by 

/dm <n> = (1 - a)-" / / o *-"dm, 

for all / G C(I <n> ). Similarly, the corresponding differential operator, H <n> — 
~~Sm~ — lx on I<n> = [0, a - rl ] , can be defined as the infinitesimal generator of the 
Dirichlet form (a<„>, V <n> ) given by 



P (ffdx = ft"a(/ofi»), for / G P<„>, 
Jo 



a<n>(f) 
where 

^<«> = {/ € L 2 (I <n> ,m <n> ) : f exists a.e. and /' G L 2 (I <n> ,dx)}. 

We define H <oa> as the operator — dm ^ - with Dirichlet boundary condi- 
tions on /<oo> = [0, oo). It is clear that the (projective system of) measures m <n> 
give rise to a measure m <00 > on I<oo> since for any / G T> <n> with support con- 
tained in [0,1], a <n> (f) = a(f) and J z fdm <n> = jjfdm. Furthermore, we 
define the corresponding Dirichlet form (a <00> , T> <00> ) by 

a<oo>(/) = supa< n> (/|/ <n> ) = lim a<n>(/|/< n >), for / G X> <00 >, 

where 

oo>) • SUp ft<n> {f\l <n > ) ^ Oo}. 



Clearly, a <00> satisfies a self-similar identity analogous to Equation (3.4 1 and its 
infinitesimal generator is H <00> . 
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3.2. The eigenvalue problem. The study of the eigenvalue problem 

(3.5) H<n>f = -*—*f = X f 

am <n> ax 

for the Sturm Liouvillc operator H <n> , equipped with Dirichlet boundary condi- 
tions on /<„>, revolves around a map p, called the renormalization map, which is 
initially denned on a space of quadratic forms associated with the self-similar set 



I <n > (or with a corresponding finite graph, see Remark 3.1 below) and then, via 



analytic continuation, on C as well as (by homogeneity) on P (C), the complex 
projective plane. (More precisely, the renormalization map is associated to the 
passage from /<„_!> to J<„> or, equivalently, from /<o> to /<i>.) As will be 



explained in §3.3, the propagator of the above differential equation (3.5| is very 
useful in producing this rational map, initially defined on C 3 , and later on, as the 
polynomial map 

(3.6) p([x, y,z]) = [x(x + S~ l y) - <T V, 5y(x + S^y) - Sz 2 , z 2 }, 

defined on the complex projective plane P 2 (C). Here, [x, y, z] denotes the homo- 
geneous coordinates of a point in P 2 (C), where (x, y, z) € C 3 is identified with 
{fix, /3y, f3z) for any j3 € C, (3 ^ 0. Clearly, in the present case, p is a homogeneous 
quadratic polynomial. 

(Note that if p is viewed as a map from P 2 (C) to itself, then one should write 
more correctly, 

p([x :y:z}) = [x(x + S~ l y) - 5~ 1 z 2 : 5y{x + S~ 1 y) - Sz 2 : z 2 }, 

where [u : v : w] — [u, v, w] denote the homogeneous coordinates of a generic point 
in the projective plane P 2 = P 2 (C). We will use this notation in §4.) 

As we shall see in §3.3, the spectrum of the fractal Sturm-Liouville operator 
H <oc> , as well as of its finite graph (or rather, here, bounded interval, see Remark 



3.1 below) approximations H <n> (n = 0, 1, 2, ...), is intimately related to the iter- 
ation of p. The spectrum of H <0> and of H <n> (n = 1,2, ...) is discrete for any 
value of a in (0, 1) (and hence, for any 7 > 4). (See, e.g., [8], Q2], 0H, [H], 
However, in the sequel, we will focus our attention on the case where a < | (or 
equivalently, 5 < 1). In that case, the spectrum of H <OD> is pure point for a < 5 
(i.e., 7 > 4), but absolutely continuous for a = | (i.e., for 7 = 4). Furthermore, 
observe that since a £ (0, 1) and 7 = a ^i_ a ) (see Equation (3.1)), we always have 



7 > 4 (and in particular, 7 > 1), independently of the above assumption according 
to which a < |. Finally, note that 7 = 4 if and only if a = ^ (i.e., 6 = 1), an 
interesting special case which will be investigated at the end of §5. 

Remark 3.1. In |44) (see also |45j and [46 ), is also provided a description of 
the renormalization map in terms of lattice or finite graph (rather than of bounded 
self-similar interval) approximations of the half-line. By necessity of concision, we 
will not discuss this matter further in this paper. The method, however, is very 
analogous to the one described in part of §4 below in the case of the infinite lattice 
SG <00> based on the bounded Sierpinski gasket SG. (In some definite sense, the 
bounded self-similar interval I = [0, 1] would play the role here of SG, while the 
half-line /<oo> = [0, 00) would be a substitute for SG <00> .) 
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3.3. The renormalization map and the spectrum of the operator. The 

renormalization map p is a function from the complex projective plane P 2 (C) to 
itself which is induced by the above homogeneous polynomial map from C 3 to itself; 



see Equation (3.6) and the discussion surrounding it. As we shall see later on, the 
spectrum of the operator H <00> defined in §3.1 is intimately connected with the 
iteration of the renormalization map. Following |44j . we next explain how the 



explicit expression for the renormalization map given in Equation (3.6) above can 
be derived by studying the propagator for the eigenvalue problem associated with 
the operator. 

We define the propagator T\(s, t) for the eigenvalue problem — dm ^ - j^f = Xf 
associated with the operator H <00> on /<oo> = [0, oo) as a time evolution function 
which, for each < s < i, is a 2 x 2 matrix with nonzero determinant such that 
the solution of the equation satisfies 



/(*) 
/'(*) 

where /' denotes the derivative of /. 



= r A (M) 



m 

/'(«) 



Using the self-similarity relations (3.2) and (3.4) satisfied by the measure m 



and the Dirichlet form a, respectively, and recalling that 7 is given by Equation 
(3.1), we obtain T <n>t \ = D a n o T~ ( n\ o D a -n for the propagator T <n>t \ associated 



with the eigenvalue problem 



d 



dm <n > dx 



f-f = Xf, where 



1 

a n 



In order to derive the expression of the renormalization map stated in Equation 
(3.6), we will consider the case when n = 1. Thus, we have r<i> \ = D a o T 7 a ° 
D a -i. Let 

o(A) 6(A)" 
c(A) d(A) ' 
We will proceed with the following calculations: 



a •= 



<1>,A 



1 

a 



o( 7 A) 6( 7 A) 
c( 7 A) d( 7 A) 





'1 " 






a- 1 





a( 7 A) a _1 6(7A)' 
ac( 7 A) rf( 7 A) 



On the other hand, we have 



t<i>,a = r A (i,a 




- 


= D s oT 


A 


oD 5 - 


1 r A 




'1 0" 




'o(A) 6(A)" 




'1 " 




o(A) 


6(A)" 




s 




c(A) d(X) 




s- 1 




c(A) 


d(A) 



a(A) 2 + ( J- 1 6(A)c(A) a(A)6(A) + ^ 1 d(A)6(A)' 
5a(X)c(X) + c(X)d(X) 5b(X)c(X) + d(X)d(X) 

Using the fact that T\ G SL2 (C) (the special complex linear group of 2 x 2 matrices) 
and hence, that a(X)d(X) — 6(A)c(A) = 1, we see that the two diagonal entries can 
be rewritten as 

'd(X)a(X) - 1 



a(A) 2 +<5- 1 6(A)c(A) = a(A) 
5b(X)c(X)+d(X)d(X) = 6d{X) 



a(X)+6- 1 
a(X)d(X) 



d(X) 



o(A) 
+5 _1 d(A) 



a(X)(a(X)+S- 1 d(X))-S-\ 
Sd(X)(a(X)+S- 1 d(X))-S. 
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We initially define the renormalization map p : C 2 — > C 2 in terms of the above 
diagonal entries as 

p(x, y) = (x(x + S~ 1 y) - 5~ x , Sy(x + 6~ 1 y) - 5) 

and the map : C — > C 2 as 

[o(A)l 



0(A) = 



d(X) 



Note that po 0(A)) = 0(7A) for all A G C. 

We now go back to the 2-dimensional complex projective space P 2 = P 2 (C), 
and note that any point [x, y, z] G P 2 is equivalent to [- , | , 1] for z ^ 0. We can 
therefore represent P 2 by 

P 2 = {(91,52,1) : («!,<&) € C 2 }U {[ar,i/,0] : (x,y) G C 2 }. 

We can then naturally define the map p, now viewed as a polynomial map from the 
complex projective plane P 2 = P 2 (C) to itself, as follows: 

(3.7) p([x, y,z}) = [x(x + 8^y) - 5~ 1 z 2 ,5y(x + S~ 1 y) - Sz 2 , z 2 }, 



which is in agreement with Equation ( |3.6[ ). 

In light of the above discussion, the invariant curve can be viewed as a map 
: C — > C 3 defined by 0(A) = (a(A), b(X), 1) and satisfying the functional equation 

(3.8) po0(A) = 0( 7 A), 

for all A G C. 



Next, we study the spectrum of the eigenvalue equation (3.5), as well as of 
its counterpart for n = oo. An attractive fixed point Xo of p is a point such that 
pxo = xq and for any other point x in some neighborhood of Xq, the sequence 
{p nx }^=o converges to xq. The basin of attraction of a fixed point is contained in 
the Fatou set of p. For S > 1, aio = [0, 1, 0] is an attractive fixed point of p. The set 

(3.9) D:={[x,y,z]:x + 6-^ = 0} 

is part of the Fatou set of p since it is contained in the basin of attraction of xq. (For 
various notions of higher-dimensional complex dynamics, we refer to the surveys 
provided in [lOj and |llj .) The set D and the invariant curve of p together 
determine the spectrum of H <n> and of H <00> . Moreover, the set of eigenvalues 
(i.e., here, the spectrum) can be described by the set 



(3.10) S := {A G C : 0(7 _1 A) G £>}, 

the 'time intersections' of the curve 0(7 _1 A) with D. It turns out that S is countably 
infinite and contained in (0, oo). We write S = {Xk}^ =1 , with Ai < A2 < ... < Xk < 
each eigenvalue being repeated accordingly to its multiplicity. (It also turns out 
that the A^'s are all simple, in this case.) Furthermore, following |28j . we call S 
the generating set for the spectrum of H <n> , with n = 0, 1, 00. 

Let S p = j p S, for each p G Z. Recall that we are assuming throughout that 
a < | (i.e., S < 1). It follows that not only the spectrum of H <n> (for n = 0, 1, 2, ...) 
is discrete but unless a = | (i.e., 5 = 1 or equivalently, 7 = 4), so is the Dirichlet 
spectrum of H <oa> . Furthermore, the spectrum of H <n> and of H <oc> can be 
deduced from that of H <0> , as will be seen in the next theorem. Finally, recall 
that we always have 7 > 4 and hence, 7 > 1. 
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Theorem 3.2 (Sabot, |47j). The spectrum of H <0> on I = /<o> is equal 
to U^Lo^P' w hHe {if a < §) the spectrum of H <00> on /<oo> — [0, oo) is equal 
t° U^l-oo Furthermore, for any n > 0, the spectrum of H <n> is equal to 

[J p —_ n Sp. Moreover, for n = 0,1,..., oo, each eigenvalue of H <n> is simple. (In 
particular, each Xj € S has multiplicity one.) 

The diagram associated with the set of eigenvalues of the operator H <QO> can 
be represented as follows: 



7 2 ^i 7 2 ^2 7 2 ^3 7 2 ^4 

7 _1 Ai 7 _1 A 2 J~ 1 X 3 7 _1 -^4 

Ai A2 A3 A4 

7A1 7A2 7A3 7^4 

7 2 Ai 7 2 A 2 7 2 A 3 7 2 A 4 



Sabot's work ( |43j - [47j ) has sparked an interest in generalizing the decimation 
method to a broader class of fractals and therefore, to the iteration of rational 
functions of several complex variables. 

For each k > 1, we denote by fk the solution of the equation H <oa> f = Xkf 
for Xk G S. In other words, fk is an eigenfunction of H <oa> associated with 
the eigenvalue Xk € S. (Note that fj~ is uniquely determined, up to a nonzero 
multiplicative constant which can be fixed by a suitable normalization.) 

Theorem 3.3 (Sabot, 07]). Assume that a < |. 
(i) Given any k > 1 and given p € Z, if fk is the normalized solution of the 
equation H <OQ> f = Xkf for Xk € S, then fk. p := fk is the solution of 
the equation H <oc> f — Xk, P f, where Afc )P := 7 p Afc and p £ Z is arbitrary, 
(ii) Moreover, if fk,p — fk° ^1 is the solution of the equation H <oc> f = Xk, p f, 
then fk,p,<n> '■= fk,p\i <n> , the restriction of fk, p to /<„>, is the solution of 
the equation H <n> f = Xk, p f- 
Finally, for each n — 0, 1, 00, {fk,p,<n> '■ k > l,p > — n} is a complete set of 
eigenf unctions of H <n> in the Hilbert space L 2 (R + , m <00> ), where R + = /<oo> = 
[0,oo). 

4. An infinite lattice based on the Sierpinski gasket 

In this section, we will show that the polynomial induced by the decimation 
method in the case of the classical (bounded) Sierpinski gasket SG can be recovered 
from the infinite lattice SG^°°^ based on the (bounded) Sierpinki gasket. 

We start with a self-similar set F = {1,2,3}, the vertices of an equilatral 
triangle, and construct an increasing sequence of finite sets F <n> by blowing-up 
the initial set F = F <0> . For instance, -F<i> is defined as the union of three copies 
of F. Namely, F = Uf =1 F<i>,i. The unbounded set F <00> , called the infinite 
Sierpinski gasket and also denoted by SG^ 00 ^ (see Figure [4 1, is the countable set 
defined as the union of all the finite sets F <n> = <E> -1 o $f o 



o 1 (F), where 



g, then the spectrum of ii<oo> is given by the closure of UJ; 
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Figure 4. An infinite Sierpinski gasket, F <co> = SG(°°\ 



(ii, %2, ■■■i n ) G {0,1,2}™, n > 0, and ($0,^15^2) are the contraction mappings, 
expressing the self-similarity of the set F: 

00 

F<oo> = [J F<n>- 
n=0 

Similarly, we can define a sequence of operators H <n> on F <n> , and the Laplace 
operator H <00> on F <oc> can be viewed as a suitable limit of the operators H <n> . 
The operators H <n> arise from the Dirichlet forms A <n> defined on £{F <n> ) = 
R F<n> , the space of real- valued functions on F <n> . To construct the Laplacian on 
the sequence F <n> , define the Dirichlet form by 

A <n> (f(x)) = J2(f(y)-f(x)) 2 , 

where / € M F< ' ,> and the sum runs over all the neighbors y of x in the finite graph 
associated with F <n> . 

The operator H <n> on L 2 (F <n> ,b <n> ) is the infinitesimal generator of the 
Dirichlet form defined by < A <n> f,g >= — J H <n> fgdb <n> , where b <n> is the 
positive measure on F <n> which gives a mass of 1 to the points in dF <n> and 2 to 
the points of F <n> \dF <n> , where dF <n> denotes the boundary of the graph F <n> . 
The sequence H <n> is uniformly bounded and we can use it to define the operator 
H<oo> on L 2 (F <QO>1 6<oo>), where b <oa> is the positive measure on F <00> defined 
as a suitable limit of the measures b <n> . 

Sabot [45 has shown that the extended decimation method which he estab- 
lished for these fractals naturally involved the complex dynamics of a renormaliza- 
tion map of several complex variables. 

Let G be the group of symmetries acting on the finite lattices F <n> , namely 
S3, which is also the natural symmetry group of the equilateral triangle and of 
F = {1,2,3}. We denote by Sym the space of symmetric (linear) operators 
on C F which are invariant under G. For every Q 6 Sym G , we can construct a 
symmetric operator Q<i> on C F<1> as a sum of the copies of Q on F < i >i j (where 
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we recall that F < i >i i, F < i >i 2 and i ? < i >i3 are the three copies of F<i>): 

3 

Q<i> = ^ Q<l>,i- 

i = l 

We define a rational map T : Sym — > Sym by 

(4.i) = (0<i>)I»f<i>. 

the trace of Q on the 'boundary' of i*<i>. In general, using the interpretation of 
symmetric matrices, it is difficult to analyze the map T, which involves the notion of 
trace on a subset. (Clearly, one can also view Q as a quadratic form on C f ; hence, 
the notation.) To avoid the complication of taking the trace of a symmetric matrix 
on a subset, the space of symmetric matrices can be embedded into a Grassmann 
algebra, in which case the map T becomes linear. 

Let E and E be two linear subspaces of with canonical basis (fj x )xeF and 
(j]x)x£F, respectively, and such that C F = E@E. We define the Grassmann algebra 
associated with this (orthogonal direct sum) decomposition of C F by 

2\F\ 

/\(£ffl£) = 0(£ffi£) Afe , 

fe=0 

where (E © E) Ak is the k-fold antisymmetric tensor product of E © E with itself, 
and \F\ — 3 is the cardinality of F . We consider the subalgebra of the Grassmann 
algebra generated by the monomials containing the same number of variables fj and 
77, namely, 

l-FI 

A=@E Ak AE Ak . 

fc=0 

We can embed Sym G into A via the injection map 



Sym G -> A 

Q h-> exp(^Qry), 

where fjQr/ = . gF QijrjiVj an( i Q : — (Qij)i,jeF> a symmetric matrix. Wc 
denote by V{A) the projective space associated with A. Let the corresponding 
canonical projection be tt : A — > V(A). The closure of the set of points of the form 
7r(exp (77(577) for Q 6 Sym , denoted by L G , is a smooth submanifold of ^(A) of 
dimension dim(Sym G ). (We refer, for example, to |42| and [?] for an introduction 
to Grassmann algebras and projective geometry.) 

With T defined by ( |4.1[ ) , we define the following homogeneous polynomial map 
of degree 3 (recall that 77^ = fj 2 =0): 

R:A -> A 
exp(fjQr]) 1 — > rfei((g<i>)| F<1> \ aF<1> exp(fjTQr,). 

This map induces a map g on the projective space L G = P 1 x P 1 such that 
g n (n(x)) = n{R n (x)) for x e tt-^L 6 ) and n > 0. (Here, P 1 =P 1 (C)=CU {00} 
is the complex projective line, also known as the Riemann sphere.) Since L G is a 
compactification of Sym G , the map g : L G — > L G is an extension of the trace map 
T : Sym G — > Sym G from Sym G to L G . 
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More precisely, note that C F can be decomposed into a sum of two irreducible 
representations: C F = Wo © W\, where W$ is the subspace of constant functions 
and Wi is its orthogonal complement. Let p\w andp|wi denote the orthogonal pro- 
jections of C F onto Wo and W\, respectively. Via the isomorphism between Sym 
and C 2 (recall that F has cardinality 3), every element Q £ Sym can be written 
as Q" '" 1 = uop\w + uip\wn where (uo,u\) £ C 2 . Then, the renormalization map 
T is defined by 

/ 3u wi ui{un+u\Y 

T{u ,ui) = ■ , — ■ 

\2u + ui oui + Uq 

This map induces another map on P 1 , denoted by g and given by 

g(i z o ■ Zi]) = N(5zi + z ) : (2z + Zx){z a + Zi)]. 

(Here, [wq : w{\ denotes the homogeneous coordinates of a generic point of the 
projective line P 1 = P 1 (C).) Indeed, upon the substitution z = the map T 
gives rise to 

g( z ) 



2« +«i 



5«i+«o 

^1(5+^) 



Ul (2^ + l)«?(aa + l) 
_ z(5 + z) 
~ (2z + l)(z + l)' 

The dynamic of the map g plays an essential role in the study of the spectrum of 
the associated symmetric operators. After having made the additional change of 
variable v = , we obtain the following equation 

3v(2v + 5) 
(3w + 3)(2i; + 3)' 

from which we recover the polynomial p(v) = v(2v + 5). This polynomial plays a 
significant role in the case of the bounded Sierpinski gasket (as shown in §2) and 
was introduced (by completely different methods) in the initial work of Rammal 
|40j and of Rammal and Toulouse |41j . We note that the dynamics of rational 
maps in higher dimension is hidden in the one-dimensional case of the Sierpinski 
gasket. 

Remark 4.1. A thorough discussion of the symplectic and super symmetric as- 
pects of the new methods developed in |45] (as well as, implicitly, in [43j and [44 ) 
to extend to the multi-variable case the classic decimation method, is provided in 
|46j . In particular, beside supersymmetry (which is translated mathematically by 
the presence of Grassmann algebras and variables) , geometric quantization and the 
associated momentum map in symplectic geometry play an important role in this 
context. 

Remark 4.2. A good review of many of the rigorously established properties 
of the Laplacian on various realizations of "the" (deterministic) Sierpinski gasket 
can be found in [53], both in the case of discrete and continuous spectra. 
( We only discuss the case of discrete spectra in the present paper; see, however, 



Remark 5.14 below.) Moreover, in [45] . the emphasis is on the study of the spectral 
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properies of random (rather than deterministic) realizations of "the" infinite (or 
unbounded) Sierpinski gasket SG^K 

5. Factorization of the spectral zeta function 

In this section, we show that the spectral zeta function of the Laplacian defined 
on a (suitable) finitely ramified self-similar set or on an infinite lattice based on this 
fractal can be written in terms of the zeta function associated with the renormaliza- 
tion map. We will focus here on the case of the Laplacian on the bounded Sierpinski 
gasket SG (as in §2 and [55j ) or on the infinite Sierpinski gasket SG^- 00 " 1 (as in §4 
and [28]), as well as on the case of fractal Sturm-Liouville differential operators 
on the half-line L <oa> = [0, oo), viewed as a blow-up of the self-similar interval 
I = [0, 1] (as in §3 and [28]). 

Definition 5.1. The spectral zeta function of a positive self-adjoint operator 
L with compact resolvent (and hence, with discrete spectrum) is given (for Re(s) 
sufficiently large) by 

oo 

(5.1) a(sH5>,r /2 , 

where the positive real numbers Kj are the eigenvalues of the operator written in 
nonincreasing order and counted according to their multiplicities. 

A. Teplyaev (|55j. see also [54]), motivated by the known identity for fractal 



strings (see Remark 5.3 below), has studied the spectral zeta function of the Lapla- 
cian on SG and, in the process, has explored interesting connections between the 
spectral zeta function and the iteration of the polynomial induced by the decimation 
method. 

Theorem 5.2 (Teplyaev, [55]). The spectral zeta function of the Laplacian on 
SG is given by 

(5.2) 

. , . . ^ > 5 ^ ( 1 3 \ A . . 5~ s / 3 1 

^W = C«, 3 (^^ 1 _3. 5 _.+ r ^J+C«, S (^^_3. 5 _ 5 - r ^ 

where R(z) := z(5 — Az) and 

(5.3) C RtZ0 (s) := lim £ (c**)"* 

zeR- n {z } 

is the polynomial zeta function of R, defined for Re(s) > 2 ^° s e 2 (where c := 5 = 
i?'(0)). Furthermore, there exists e > such that Ca m ( s ) has a meromorphic con- 
tinuation for Re(s) > —e, with poles contained in ^ f^f , los : Ti G Zj. 

Remark 5.3. Theorem \5.S\ extends to the present setting of analysis on certain 
fractals |25j (or 'drums with fractals membrane', see, e.g., [7], [9], |13j — |lTj . [20] . 
[26] -[28], [31] -[32], [35], 00] -05], [50], [51], [53] -[55]) and provides a dynamical 
interpretation of the factorization formula obtained by the second author in [30 , 
13 1 j for the spectral zeta function C, v (s) of the Dirichlet Laplacian associated with a 
fractal string (a one- dimensional drum with fractal boundary [29 -{36 ): 

(5.4) U*) = CM • C», 
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where £(s) is the classic Riemann zeta function and Cg( s ) * s (the meromorphic con- 
tinuation of) the geometric zeta function of the fractal string (whose poles are called 
the complex dimensions of the string and help describe the oscillations intrinsic to 
the geometry and the spectrum of the string, see [36 ). See also ([36], §1.3) for 



a discussion of Equation (5.4 1 and, for example, [301 . |31j . |33j . [34j . along with 
much of |36) . for various applications of this factorization formula. Finally, we 
note that in [55], £(s) is reinterpreted as the polynomial zeta function of a certain 
quadratic polynomial of a single complex variable and hence, the factorization for- 



mula (5.4) can also be given a complex dynamical interpretation (in terms of the 



iteration of the renormalization map). This issue was revisited in 28j, and, as 



a result, formula (5.4) can also be interpreted in terms of multi-variable complex 



dynamics. (See [28 , along with Theorem 5.12 below 



Remark 5.4. In [9], shortly after the completion of [54], [55], G. Derfel, P. 
Grabner and F. Vogl, working independently on this question and motivated in part 
by the results and conjectures of |26j and |31j , have obtained another interpretation 



of the geometric factor of the factorization formulas (5.2) and (5.4), expressed in 
terms of the multiplicities of the eigenvalues. In the process, they have shown that 
the spectral zeta function £a (s) has a meromorphic continuation to all ofC More- 
over, they have proved further cases of a conjecture of |26j . [31] , according to which 
the asymptotic second term in the spectral counting function of the Laplacian on the 
Sierpinski gasket SG and other lattice self-similar fractals is truly oscillatory (or 
equivalently, the corresponding periodic function obtained in |26) is not constant). 
[Compare with analogous conjectures and results stated or obtained in |29j - [36] for 
drums with fractal boundary (instead of drums with fractal membrane).] 



Later on, in [28], we extended this result about the factorization of the spectral 
function of the Laplacian on the finite (or bounded) gasket SG to the infinite (or 
unbounded) Sierpinski gasket SG^°°^ and to the renormalization maps of several 
complex variables associated with fractal Sturm-Liouville operators. 

The Dirac delta hyperfunction on the unit circle T is defined as St = [8j,Sj] = 
[jzr, —zr]- It consists of two analytic functions, Sj : E — > C and Sj : C\E — > C, 
where E — {z G C : \z\ < 1 + ^} for a large natural number N. In other words, 
a hyperfunction on T can be viewed as a suitable pair of holomorphic functions, 
one on the unit disk \z\ < 1, and one on its exterior, \z\ > 1. (See, for example, 
( |19j . §1.3) and f |37| . §3.3.2) for a discussion of various changes of variables in a 
hyperfunction. See also those two books |19j . |37j . along with |48j . |49j . the original 
articles by M. Sato, for an overview of the theory of hyperfunctions. Moreover, see 
|52j for a detailed discussion of St and, more generally, of hyperfunctions on the 
unit circle T.) 

Theorem 5.5 ( 28] , Lai and Lapidus). The spectral zeta function Ca<°°) of the 
Laplacian A^ 00 ) on the infinite Sierpinski gasket SG^ 00 ^ is given by 

(5.5) Ca(-)(«)=Ca m (*)*t(5-5), 

where St is the Dirac hyperfunction on the unit circle T and Ca is the spectral 
zeta function of the Laplacian on the finite (i.e., bounded) Sierpinski gasket SG, as 



given and factorized explicitly in Equation (5.2) of Theorem 5.2 
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In the case of the Sturm-Liouville operator on the half-line, we introduce a 
multi- variable analog of the polynomial zeta function occurring in Equation (5.3) 
of Theorem 15.21 

We define the zeta function of the renormalization map 



Definition 5.6 ( 
p to be 



(5.6) 



Us) 



OO 

E E 

p=o {xec-. pp(0(7-(p+ 1 U))ei5} 



(7 p A)" 



for Re(s) sufficiently large. 

Remark 5.7. Recall from §3 that in the present situation of fractal Sturm- 
Liouville operators, the renormalization map p is given by Equation (3.6) and that 
the 'renormalization constant' (or 'scaling factor') 7 is given by Equation ( |3.1[ ). 
(Also, see the functional equation (3.8) defining the invariant curve <f>, as well 
as Equations (3.9) and (3.10) defining D and the generating set S, respectively.) 
Furthermore, recall from the end of §3.2 that we assume that a < | (i.e., 5 < 1) 
in order to ensure the discreteness of all the spectra involved^ and that we always 
have 7 > 4; in particular, 7 > 1, and 7 = 4 if and only if a = ~. 



We consider the factorization formulas associated with the spectral zeta func- 
tions of the sequence of operators H <n> = — dm ^ - , starting with H <0> on [0,1], 
which converges to the Sturm-Liouville operator H <co> on [0, 00). In light of Def- 
inition [5T] and Theorem 3.2 given any integer n > 0, the spectral zeta Ch<„> ( s ) of 
] is initially given by ( H<n> (s) = J2\ e s E^L-„(7 P A)~5 , for Re(s) 



H 



<n> 



on [0,c 
large enough. 



Theorem 5.8 (|28J. The zeta function ( p (s) of the renormalization map p is 
equal to the spectral zeta function of H <0> , 

00 

Ch<o>(^EE(7 P A)^, 

\£Sp=0 

or its meromorphic continuation thereof: Cp( s ) = Ch <0> ( s ) ■ (An expression for 
Ch <0 > (s) is given by the n — case of Proposition \57^ just below; see Equation 
§M of Remark\5A(\) 



(5.7) 



&<«>(«) 



<s(s), 



Proposition 5.9 ([28 ). For n>0 and Re(s) sufficiently large, we have 

(7") f 
1-7" 

where Cs(s) is the geometric zeta function of the generating set S. Namely, £s(s) :— 
Y^jLi(^j) _5 (for Re(s) large enough) or is given by its meromorphic continuation 
thereof. 

Remark 5.10. In particular, in light of Theorem \5.£\ we deduce from the n = 
case of Proposition \5.9\ that 

(5-8) Us) = Ch <0> (s) = ^^jCs(s). 



1 



7 



Except when a - 



|, in which case the spectrum of H< oc> is continuous. 
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In the case of the operator H <oc> , the asymptotic behavior of the spectrum 
led us naturally to using the notion of delta hyperfunction. Indeed, according to 
Theorem|3.2| a part of the spectrum of H <oa> tends to while another part tends 



to oo. If one mechanically applies Definition 5.1 one then deduces that Cff<oo>( s ) 
is identically equal to zero, which is clearly meaningless. Indeed, the geometric 
factor in the factorization of Ch <00> (s) is equal to the sum of two geometric series 
(converging for Re(s) > and for Re(s) < 0, respectively), and this sum is itself 
identically equal to zero (except for the fact that one of the two terms in the sum 
is not well defined, no matter which value of s G C one considers). Fortunately, 
there is a satisfactory resolution to this apparent paradox. More specifically, we 
have discovered that the geometric part of the product formula for the spectral zeta 
function Ch <00> carL be expressed in terms of the Dirac delta hyperfunction Sj on 
the unit circle. 

Theorem 5.11 ([28]). Assume that a < h. The spectral zeta function Ch <00> 
is factorized as follows: 

(5-9) = Cs(s) ■ S T (j-i) = ( p (s)(l - 7~ 4 )M7~ 4 )- 



Furthermore, we have shown in 28j that the zeta function associated with the 
renormalization map coincides with the Riemann zeta function £(s) for a special 
value of a. When a = \ (or equivalently, 5 = 1 and so 7 = 4), the self-similar mea- 
sure m coincides with Lebesgue measure on [0, 1] and hence, the 'free Hamiltonian' 
H = H <0> coincides with the usual Dirichlet Laplacian on the unit interval [0, 1]. 



Recall that Cp( s ) = Cfr< > ( s ) 0°Y Theorem 5.8 1 and that the spectrum of H 



l <o> 



(i.e., of the Dirichlet Laplacian on [0, 1]) is discrete and given by 7r 2 j 2 , for j = 1, 2, .... 
Hence, Cp( s ) = Ch <0> ( s ) = t ~ s (,(s), where £(s) is the Riemann zeta function. 

Theorem 5.12 ( |28j ). When a = g, the Riemann zeta function £ is equal (up 
to a trivial factor) to the zeta function £ p associated with the renormalization map 
p <mP 2 (C). More specifically, we have 

(5.10) C(«) = ^Cp(s) = 7-^37 Cs(s), 



where £ p is given by Definition 5.6 (or its mermorphic continuation thereof) ano 



the polynomial map p : P 2 (C) -> TF(C) is given by Equation (|3.6[) with a = i (and 



hence, in light of (3.1), with 8 = 1 and 7 = 4): 

(5.11) p([x, y, z]) = [x(x + y)- z 2 , y(x + y)- z 2 ,z 2 }. 

This is an extension to several complex variables of A. Teplyaev's result 55J, 
which states that the Riemann zeta function can be described in terms of the zeta 
function of a quadratic polynomial of one complex variable (as defined by Equation 



(5.3)) 



Remark 5.13. An interesting open problem is to determine what happens if 
we consider Ch <<X3> (s) instead of Ch <0> ( s )> s ^ when a = |. In that case, the 
spectrum of H <oa> is continuous and a suitable interpretation has to be found for 
Ch <00 >( s )> i n terms of a properly defined spectral density of states. We leave the 



investigation of this open problem for a future work. (See also Remark 5.14 just 
below.) 
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Remark 5.14. It would be interesting to obtain similar factorization formu- 
las for more general classes of self-similar fractals and fractal lattices, both in the 
present case of purely discrete spectra or in the mathematically even more chal- 
lenging case of purely continuous (or mixed continuous and discrete) spectra. The 
latter situation will require an appropriate use of the notion of density of states, of 
frequent use in condensed matter physics (see, e.g., PQ, [38], [39]) and briefly dis- 
cussed or used mathematically in various related settings (involving either discrete 
or continuous spectra) in, e.g., [17j . [20] , [26 and especially, |43] , [44] . |45j . 
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